Characterization of the Emergence of Order in an Oscillated Granular Layer 
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The formation of textured patterns has been predicted to occur in two stages. The first is an 
early time, domain-forming stage with dynamics characterized by a disorder function 5{I3) ~ t'"^ , 
f*"^ I with as = |/9; this decay is universal. Coarsening of domains occurs in the second stage, in which 

. (5(/3) ~ t~''^ , where tJi is a nonlinear function of P whose form is system and model dependent. Our 

' experiments on a vertically oscillated granular layer are in accord with theory, yielding ctb « 0.5/3, 

04 , and (tl a nonlinear function of /3. 

PACS number(s); 05.70.Ln, 82.40. Ck, 47.54.+r 

The formation of textured spatial patterns has been studied in laboratory experiments |l| and model systems [^j . 
QQ ' Several quantities have been used to describe the development of patterns represented by a scalar field v{k) with a 
typical wavevector feg, including the structure factor [^|-|7j, density of topological defects and a recently introduced 
family of characterizations called the disorder function, 

.5 ■■ ' 1^'- ^r(i"(x)i>^ • 

where (|z)(x)|) denotes the mean of |f(x)|, and 5(/3) ( < /3 < 2) has been normalized to be scale invariant 
d{P) describes configuration-independent aspects of textures and their formation, aspects that are independent under 
repetition of the experiment. It can be used to study multiple aspects of patterns just as generalized dimensions pO| ] 
and singularity spectra can be used to describe multiple aspects of strange attractors. In experiments presented 
in this paper, we characterize the time evolution of patterns in a vertically oscillated granular layer using S{P). 
' Before describing our experiments, we review numerical and analytical work on the coarsening of patterns after 
lO ! a quench from an initial featureless (or noisy) state. Most studies have focused on solutions u{r,t) to the Swift- 
Hohenberg equation [||, 

^ = {e-{A + kir)u-u^-,.{Vur+rj{r,t), (2) 

(—| ^ where u(r, t) is a two-dimensional scalar field, e is the distance from pattern onset, v is the strength of a non-variational 
term and 77 a random forcing term such that {r]{r, t)ri{r' ,t')) — 2F6{r — v')6{t — t'), where F controls the strength 
of the noise. 

For the formation of textures in (^), the width of the structure factor S{t) (i.e., the width of the peak in the 
azimuthal average of (u(k, t)-u(— k, t))), has been shown to decay in two distinct stages |Q. S{t) ~ is obeyed until 
the peak amplitude of the field u(x, t) saturates, beyond which time the pattern coarsens and the decay becomes 
slower. For e — 0.25 and 1^ = 0, Cross and Meiron H, Elder et al. jH, and Hou et al. M found that in this second 
region S{t) decreased as t~5 when F — 0, and as when F ^ Schober et al. found that for F = and 

z/ = 0, S{t) ~ ; the discrepancy with earlier results could be due to the one-dimensionality of their model. 

The characterization of pattern formation in model systems using S{f3) has also shown the presence of two stages: 
the emergence of domains characterized by 6(0) ~ with asiP) — hp, and a slower coarsening behavior p2[ . 

The relaxation exponent during coarsening depends on the value of z/ in (g) , and is thus expected to be system and 
model dependent [p^ . We report analogous behavior in the present laboratory study of pattern formation in an 
oscillated granular layer and present additional differences between the two stages. 

Our experiments generate patterns in a layer of 0.165 mm bronze spheres contained in a vertically oscillated 
circular container with diameter 14 cm |lj]. The layer is four particle diameters deep, and the cell is evacuated to 4 
Pa to avoid any hydrodynamic interaction between the grains and surrounding gas. The control parameters are the 
frequency / of the sinusoidal oscillations and the peak acceleration of the container, F = (27r/)^A^/(7, where A is the 
maximum amplitude of the oscillation and g is the gravitational acceleration. As / and F are varied, a variety of 
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textures including striped, square or hexagonal planforms are observed |TJ]. Our analysis reported here is restricted 
to patterns with square planforms. In the region of the phase diagram studied, square patterns appear for increasing 
control parameter at F « 2.75; the bifurcation is subcritical. The geometry of the circular container allows relaxation 
to an almost perfect square array through wavelength adjustment of the pattern at the container wall over a distance 
of less than one wavelength . 

The granular surface is illuminated with a ring of LEDs surrounding the cell. The light is incident at low angles and 
the scattering intensity is a nonlinear function of the height of the layer; scattering from peaks (valleys) creates bright 
(dark) regions. The images are collected at the driving frequency and the acceleration of the container is monitored 
during each run. F is suddenly increased from its initial value of 2.2, where no discernible structure is observed. As 
the grains are not in contact with the container for part of their motion, we assume that the initiation of the quench 
occurs at the first layer-plate collision after the change of control parameter. The uncertainty in the time origin is 
the dominant source of error in our measurements. 

The top row of Fig. 1 shows that local square domains emerge and coarsen to a final, almost perfect, square array. 
The bottom row shows this process in Fourier space. A repetition of the experiment would lead to similar, but not 
identical, intermediate states. Our aim is to study configuration-independent aspects of this relaxation, and to analyze 
their dependence on the control parameters / and F. 

Patterns such as those shown in the top row of Fig. 1 can be represented by a discrete sampling of a smooth 
scalar field w(x). The values of the field are known on a (typically square) grid, but the analytical form of the field is 
unknown. The ingredients used to deduce the form of the disorder function in (|l|) are its invariance under arbitrary 
rigid motions of the texture and the nature of the local planform. Local deviations of a pattern from squares (due 
to curvature of the contour lines [^) contribute to S{f3) through the Laplacian, while variations of the size of squares 
contribute via the choice of a "global" fcg, which is obtained from the field ^(x) by minimizing the value of ^(1) [ p^ . 
Unlike the information contained in the structure factor, |(A -I- fcg)w(x)| is a local density of irregularities in the 
texture, and hence distinct "moments" /3 can be used to quantify multiple aspects of the disorder density. 

The images shown in Fig. 1 have sharp changes at the edges which lead to high frequency contributions in their 
Fourier spectra. Their removal through simple filtering causes contamination of the pattern near the edges and leads 
to error in calculating S{l3). We use a method of noise filtering that involves extending the image to a periodic one 
using "Distributed Approximating Functionals" (DAFs) jl^JT^. Fourier filtering can then be used on the extended 
image to eliminate high frequency noise and undesirable harmonics [ p^ . The evaluation of S{P) requires an accurate 
estimation of Aw(x), which typically amplifies any noise present in the discrete, digital experimental data. A method 
for this calculation has been presented in 

The behavior of 6{1) for the relaxation of Fig. 1 is shown by the symbols o in Fig. 2. The initial formation of 
local rectangular domains and the final coarsening correspond to distinct power law decays of (5(1). The transition 
coincides with the saturation of the peak amplitude; i.e., nonlinear effects are negligible during domain formation 
(typically 4-5 oscillations) and become relevant during coarsening 0. 

During the initial stage of pattern formation S{1) ~ ^-o.49±o.02^ This dynamical scaling is similar to that describing 
the deca y o f the width of the structure factor [Q and is related to the rate of domain growth in phase ordering 
kinetics po| , ^ . 

Since nonlinear effects are negligible during domain formation, the evolution can be modeled by (^ with the removal 
of the nonlinear and stochastic terms. Numerical integration starting from states consisting of random or Gaussian 
noise shows that / (Px\u{r,t)\ ^ e'^^t~i and ^ (Px\{/\ + k'^)u{Y,t)\ ~ e'^*t~3; consequently (5(1) ~ t~2. It has been 
shown analytically that the width of the structure factor for evolution of u(r, t) decays like S{t) ^ t^i until the peak 
amplitude saturates providing further evidence for the interpretation that the spatiotemporal dynamics is linear 
during the first stage. 

Furthermore, during the domain forming stage, moments of the disorder function decay as S{f3) ^ t^'^E(P) ^ where 
GEiP) ~ \P (see Fig. 3(a)). Analogous behavior can also be seen by numerical integration of the linear terms in (|^). 
The linearity of (t_e(/3) and the value of the proportionality constant suggest that multiple aspects of textures, such as 
structure factor, curvature of contours and defect densities decay via a single mechanism during the domain forming 
stage. 

The initiation of domain coarsening coincides with the saturation of peak heights of the granular layer (see Fig. 
2); i.e., the latter stages of pattern formation correspond to nonlinear spatiotemporal dynamics of the field 0]. The 
observed scaling of the disorder function is more complex. For the evolution shown in Fig. 2 at F = 2.8, 5{\) ~ t"^-^^ 
p2| . This is close, but not identical to the rate of relaxation of the structure factor . Even though the moments 

5{(3) decay (approximately) as power laws t^'^'^^^\ the exponent is not linearly related to (3 as during the early phase, 
see Fig. 3(b). There is a tendency towards slower relaxation in Fig. 3(b) (compared to Fig. 3(a)) for larger values 
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of (3. Since large values of (3 preferentially weight pattern defects, this is consistent with the slower relaxation of the 
density of defects js) than that of the structure factor |§,||j7) . 

The nonlinearity of (7l(/3) implies that pattern relaxation occurs on more than one length scale. If the "envelope" 
of the texture depends on a single length scale L{t), the field can be locally expanded as m(x, t) — Re [e''^"A(X, t)\ , 
where X = x/L{t); thus 

{A + ki)u{x,t) = Re^e'^'-^i^k ■ y^A - ^AxA) . 

Since L(t) ^ during domain coarsening, the last term can be neglected, leading to 

|(A + fco>(x,i)|«^k.V^~ J^. 

Consequently S{P) ^ 1/ L^{t) ~ ((5(1))'^. Thus the nonlinearity of aL{P) implies that the relaxation during coarsening 
occurs on multiple time scales. 

Next, we briefly consider changes in the behavior of the disorder function as the experimental system is driven 
further away from the onset of patterns. The decay of S{P) during domain formation remains unchanged, but the 
decay rate in the second region decreases with increasing F. Similar behavior has been observed with the addition of 
i.in§0]. 

We have shown that the formation of texture in a vertically oscillated granular layer occurs in two distinct stages. 
During the initial stage the spatiotemporal dynamics is essentially linear and the disorder function obeys a universal 
power law S{0) ~ f-'^^iP) ^ with <je{P) — \P; this simple behavior is also observed in the linearized Swift- Hohenberg 
equation. Nonlinearity of spatiotemporal dynamics becomes relevant during domain coarsening and ^ t~"^''^'^ 
where ctl(/3) is a nonlinear function of j3. The exponent (tl{I3) is model and parameter dependent p2[ . Such non- 
universal, configuration independent characteristics of pattern formation can be used to determine the validity and 
limitations of model systems [|3| . 

Although dynamical scaling has been reported in the formation of patterns in model systems |4||5|,|8|,|T^ , ours is 
the first reported observation of stages exhibiting trivial (i.e., a single scaling index such that S{f3) ~ (6{1))^) and 
nontrivial scaling during pattern formation in an experimental system. The methods introduced here are expected to 
have applications in studying other aspects of textures, such as the quantitative description of patterns in magnetic 
bubble material p^ , p5t . 
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(DJK). 
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FIG. 1. Snapshots showing the emergence of a square spatial pattern in a granular layer at / = 27 Hz and T = 2.8; the 
times given in the upper left corner of each image are in units of container oscillation periods. Each image in the first row is of 
the central 8 cm of the 14 cm diameter circular container. Each image in the second row is the Fourier transform of the image 
above it. In the first row, the first three frames show the emergence of local domains from a uniform background, and the last 
three show the slower coarsening of these domains to an almost perfect square array. 
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FIG. 2. The time evolution of tlie disorder function 5(1) for square patterns at three different final container accelerations, 
r = 2.8 (o), r = 3.0 (O), and F = 3.2 (+). Also shown is growth of the pattern amplitude at V = 2.8 (•). Rapid early growth 
in the domain-forming phase is followed by a saturation of the amplitude in the coaxsening phase. Each curve is an average of 
10 runs at the same control parameters. The error bar at late times for o shows typical variation between distinct runs. The 
error bars at early times are the size of the symbols. 



FIG. 3. The slopes aE{P), ffL{f3) of the curves 5{f3) during (a) domain formation and (b) coarsening. The results are for a 
single run at control parameters f — 27 Hz and F = 3.0. The bold straight lines are drawn to guide the eye. In (a), the bold 
line has slope = 1/2 and in (b), the bold line has slope = 0.12. 
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